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ABSTRACT 

The  spectral  energy  balance  of  waves  in  the  surf  zone 
is  examined  with  extensive  measurements  from  the  Duck94 
experiment.  Cross-shore  energy  flux  gradients  are 
estimated  from  spectra  observed  with  closely  spaced 
pressure  sensors.  Nonlinear  energy  exchanges  between 
different  wave  components  in  the  spectrum  are  estimated 
from  observed  bispectra  based  on  Boussinesq  theory  for 
near-resonant  triad  interactions.  Dissipation  of  wave 
energy  in  the  poorly  understood  breaking  process  is 
inferred  as  the  residual  term  in  the  spectral  energy 
balance . 

Analysis  of  the  spectral  energy  balance  shows  that 
large  decreases  in  energy  flux  observed  at  the  dominant 
wave  frequencies  as  waves  break  over  a  sand  bar  are 
closely  balanced  by  nonlinear  energy  transfers  to  higher 
frequencies.  That  is,  the  decay  of  the  spectral  peak 
within  the  surf  zone  is  a  result  of  nonlinear  energy 
transfers  rather  than  direct  dissipation.  At  higher 
frequencies,  observed  energy  flux  gradients  are  small  and 
do  not  balance  the  nonlinear  transfers  of  energy  to  high 
frequency  components  of  the  spectrum.  This  analysis 
suggests  that  the  spectrum  is  saturated  at  high 
frequencies,  and  thus,  the  energy  that  cascades  through 
nonlinear  interactions  to  higher  frequencies  is 
dissipated  in  the  high-frequency  tail  of  the  spectrum. 
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I.    INTRODUCTION 

The  lack  of  accurate  models  for  the  transformation  of  waves  across  the  surf  zone 
is  the  greatest  deficiency  is  in  our  ability  to  predict  surf  characteristics  (e.g.,  spilling  of 
plunging  breakers),  near-shore  circulation  (e.g.,  longshore  currents,  undertow,  rip 
currents),  and  sediment  transport  (e.g.,  beach  erosion,  accretion  and  the  formation  of  sand 
bars).  The  properties  of  ocean  surface  waves  on  beaches  are  affected  by  changes  in  water 
depth,  the  bottom  boundary  layer,  strong  nonlinear  dynamics  in  shallow  water,  and  surf 
zone  wave  breaking.  As  the  water  depth  decreases,  waves  refract  towards  the  beach.  As 
a  result,  propagation  directions  in  the  surf  zone  are  typically  close  to  normal  incidence, 
and  the  shoaling  evolution  of  waves  can  be  reasonably  well  described  by  1 -dimensional 
models  for  uni-directional  waves. 

Turbulence  associated  with  bottom  friction  causes  a  gradual  decay  of  waves  in 
shallow  water,  that  is  believed  to  be  important  to  the  evolution  of  wave  spectra  across 
shallow  continental  shelves.  On  beaches  that  typically  span  only  O(10)  wavelengths, 
bottom  damping  is  usually  neglected.  The  dominant  physical  processes  that  affect  wave 
transformation  on  beaches  are  believed  to  be  nonlinear  interactions  and  wave  breaking. 
Nonlinear  wave-wave  interactions,  weak  in  the  open  ocean,  are  strongly  enhanced  on 
beaches  owing  to  the  weakly  dispersive  nature  of  surface  gravity  waves  in  shallow  water. 
Energy  is  continuously  exchanged  in  near-resonant  interactions  between  all  triads  of  wave 
components  with  frequencies  fl5f2,f3,  that  obey  the  interaction  rule  f x  ±  f2  ±  f3  =  0  (Freilich 
and  Guza,  1984).  These  so-called  triad  interactions  not  only  cause  a  broadening  of  the 
wave  spectrum  (i.e.,  energy  transfers  from  the  dominant  waves  to  higher-and  lower 
frequencies),  but  are  also  responsible  for  the  dramatic  transformation  of  wave  shapes  from 
nearly  symmetric  sinusoidal  deep  water  swells  to  the  pitched-forward  wave  shape  of  near 
breaking  waves  on  beaches  (Elgar  and  Guza,  1985a). 

Models  for  wave  shoaling  transformation  outside  the  surf  zone  are  usually  based 
on  nondissipative  Boussinesq  equations  for  weakly  nonlinear,  weakly  dispersive  waves 
(Peregrine,  1967).  These  models  include  time-domain  formulations  (Madsen  et  ai,  1991), 
discrete  frequency  domain  formulations  (Freilich  and  Guza,  1984;  Liu  et  ai,  1985)  and 
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stochastic  formulations  that  predict  the  evolution  of  continuous  wave  spectra  (Herbers  and 
Burton,  1997;  Norheim  et  al.,  1998).  Extensive  comparisons  to  field  data  confirm  that 
Boussinesq  models  predict  accurately  the  shoaling  evolutions  of  waves  outside  the  surf 
zone  (Freilich  and  Guza,  1984;  Elgar  and  Guza  1985b;  Norheim  et  al.,  1998). 

Early  surf  zone  models  [Collins,  1970;  Battjes,  1972;  Goda  1975]  simply 
parameterized  the  heights  of  breaking  waves  in  the  terms  of  the  local  water  depth.  Battjes 
and  Janssen  [1978]  introduced  a  frequency-integrated  energy  balance  model  with  the 
dissipation  sink  term  based  on  the  energy  loss  in  a  turbulent  bore  and  a  simple  breaking 
criterium.  Thornton  and  Guza  [1983]  refined  this  model  by  incorporating  realistic  wave 
height  statistics.  Extensive  comparisons  with  laboratory  and  field  data  for  a  range  of 
beach  profiles  and  incident  wave  conditions  shows  that  the  energy  balance  models  of 
Battjes  and  Janssen  and  Thornton  and  Guza  predict  accurately  the  variations  of  wave 
heights  in  the  surf  zone  (Battjes  and  Janssen,  1978;  Thornton  and  Guza,  1983;  Battjes 
and  Stive,  1985).  Although  wave  height  decay  in  the  surf  zone  is  predicted  accurately  by 
these  bulk  models  for  total  wave  energy,  the  associated  dramatic  changes  in  the  wave 
spectrum  resulting  from  nonlinear  interactions  are  not  described  by  these  models,  and 
therefore  the  characteristics  of  breaking  waves  (e.g.,  spilling  vs  plunging  breakers;  narrow 
vs  broad  spectra)  are  not  predicted. 

Recently,  Boussinesq  models,  that  rigorously  account  for  nonlinear  interactions, 
were  extended  into  the  surf  zone  by  including  a  heuristic  dissipation  term  in  the  shoaling 
evolution  equations  (Mase  and  Kirby,  1992;  Kaihatu  and  Kirby,  1995;  Eldeberky  and 
Battjes  1995;  Chen  et  al.,  1997).  As  virtually  no  information  exists  on  the  spectral 
characteristics  of  dissipation  rates  in  the  surf  zone,  simple  parameterizations  were  used. 
Mase  and  Kirby,  (1992)  and  Kaihatu  and  Kirby,  (1995)  used  a  frequency  dependent 
dissipation  function  that  combines  a  constant  term  that  reduces  spectral  levels  uniformly 
in  proportion  to  the  energy  at  each  frequency  with  a  frequency-weighted  (f2)  term  that 
dissipates  the  higher  frequency  components,  using  an  adjustable  coefficient  to  control  the 
relative  importance  of  the  two  terms.  Eldeberky  and  Battjes  (1995)  used  a  frequency 
independent  formulation,  in  which  dissipation  is  proportional  to  the  energy  at  each 


frequency.  Chen  et  ai,  [1997]  showed  that  the  evolution  of  wave  shapes  across  the  surf 
zone  is  predicted  more  accurately  by  models  with  a  dissipation  term  that  increases  with 
frequency  than  models  using  a  frequency  independent  dissipation  term. 

To  obtain  a  better  understanding  of  the  spectral  characteristics  of  energy 
dissipation  within  the  surf  zone,  the  spectral  energy  balance  of  shoaling  and  breaking 
waves  is  examined  in  this  study  using  extensive  field  observations  from  a  natural  sandy 
ocean  beach.  Neglecting  alongshore  depth  variations,  directional  spreading  of  waves  and 
wave  reflection  from  shore,  the  energy  balance  can  be  expressed  as 


F<f>-WsJf>  o> 


where  the  left-hand  side  is  the  cross-shore  gradient  of  the  energy  flux  F(f).  The  source 
term  Snjj)  represents  the  net  nonlinear  transfer  of  energy  to  waves  with  frequency  /  owing 
to  triad  interactions.  The  other  source  term  Sds(f)  represents  energy  dissipation  during 
wave  breaking.  The  energy  flux  gradients  are  evaluated  from  closely  spaced  instruments. 
The  energy  transfers  in  near-resonant  triad  interactions  are  a  function  of  the 
phase-relationships  between  the  three  interacting  components,  and  thus  depend  on  higher- 
order  statistics  of  the  wave  field.  A  stochastic  formulation  of  the  Boussinesq  equations 
yields  the  following  expression  for  Sn(f)  (Herbers  and  Burton,\991\  Norheim  et  ai,  1998) 


S  (f)=IM{b(f))  (2a) 

nl 
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where  IM{  }  indicates  the  imaginary  part.  The  third-order  bispectrnm  B(f,,f2)  describes 
the  degree  of  coupling  and  the  phase  relationship  in  triads  of  wave  components  with 
frequency  fl9f2  and  fj+f2  {Hasselmann  et.  al,  1963).  In  shallow  water,  the  bispectrum 
evolves  strongly  and  describes  in  a  statistical  sense  the  shapes  of  shoaling  waves  (e.g., 
Elgar  and  Guza,  1985b;  Elgar  et  al,  1990a).  The  nonlinear  source  term  Snl(f)  in  the 
energy  balance  involves  an  integral  b(f)  of  Bifj,/^  over  all  triads  that  include  a  component 
with  frequency  f. 

In  the  present  study  the  terms  in  the  energy  balance  equation  (1)  were  examined 
with  measurements  from  a  cross-shore  array  of  pressure  sensors  deployed  near  Duck  N.C., 
in  1994.  The  energy  flux  gradient  FJf)  and  nonlinear  source  term  Snl(f)  were  estimated 
from  measured  spectra  and  bispectra  at  14  cross-shore  positions.  The  associated 
dissipation  rate  SJf)  was  then  inferred  as  the  residual  term  in  equation  1 .  The  field  data 
and  analysis  are  described  in  section  2. 

Example  observations  of  the  energy  balance  in  a  range  of  conditions  are  presented 
in  section  3,  followed  by  a  statistical  analysis  of  the  entire  2  months  observations  in 
section  4.  In  the  surf  zone,  large  negative  values  of  Fx(f)  in  the  energetic  part  of  the 
spectrum  are  shown  to  be  approximately  balanced  by  Snl(f).  These  results  show  that  the 
energy  of  the  dominant  waves  is  not  directly  dissipated  in  the  surf  zone,  but  transferred 
to  higher  frequencies  through  nonlinear  interactions  where  the  energy  is  presumably 
dissipated.  The  results  are  summarized  in  section  5. 


II.  FIELD  DATA  AND  ANALYSIS 

Detailed  field  measurements  of  the  evolution  of  shoaling  and  breaking  waves  on 
a  natural  ocean  beach  were  acquired  during  the  Duck  94  experiment  at  the  U.S.  Army 
Corps  of  Engineer  Field  Research  Facility  located  near  Duck,  N.C.,  on  a  relatively  straight 
barrier  island  exposed  to  the  open  Atlantic  ocean.  A  dense  cross  shore  transect  of  14 
pressure  sensors  was  deployed,  extending  from  the  shore  line  to  about  5  m  depth  (350  m 
from  shore) (Figure  1).  Collocated  with  each  pressure  sensor  was  a  sonar  altimeter  that 
provided  accurate  estimates  of  the  location  of  the  sea  floor.  The  sample  frequency  of  all 
sensors  was  2  Hz.  High  quality  data  were  collected  nearly  continuously  during  the 
months  of  September  and  October  1994  (Elgar  et  ai,  1997). 

The  closely  spaced  pressure  sensors  allowed  for  direct  estimates  of  the  advection 
term  FJJ)  and  the  nonlinear  source  term  S„,(f)  in  the  energy  balance  equation  (1).  A 
surface  elevation  spectrum  E(f)  and  bispectrum  B(fj,f2)  with  resolution  of  .0078  Hz  were 
estimated  from  each  hour-long  pressure  record  using  the  shallow  water  approximation. 

The  spectral  energy  balance  was  evaluated  at  the  midway  point  between  adjacent 
instruments  using  a  simple  finite  difference  approximation.  The  energy  flux  gradient  Fx(f) 
was  approximated  by  (using  the  shallow  water  approximation  of  the  group  speed): 


E(f){ghf-E(f){gh)2 

F(/)«— - —  (3) 

x  x  -x 

2       1 


where  xlrx2(>X!)  are  the  cross-shore  positions  of  a  pair  of  adjacent  sensors,  hj  and  h2  are 
the  corresponding  water  depths,  Ej(f),  and  E2(f)  are  the  corresponding  spectra. 

The   nonlinear   energy  transfer  term   (Snl)   at  the   same  midway   point   was 
approximated  by: 

SJfl^ibfl+bJU))  (4) 


with  b,(f)  and  b2(f)  the  bispectral  integral  b(f)(Eq.2b)  evaluated  at  the  instrument  locations 
xx  and  x2. 

Significant  weather  changes  throughout  the  Duck  94  experiment  provided  ample 
sampling  of  a  multitude  of  surf  conditions  including  two  major  Noreaster  storms  with 
maximum  significant  wave  heights  of  2.5  m  and  3.8  m  (Figure  2)  and  surf  zones 
extending  across  the  entire  instrumented  transect.  The  data  set  also  includes  calm  periods 
with  significant  wave  heights  as  small  as  .2  m  and  virtually  no  breaking  across  the 
transect.  Observed  wave  spectra  include  narrow  spectra  of  long  period,  remotely 
generated  swells,  broad  spectra  of  locally  generated  seas  and  mixed  swell-sea  spectra. 
The  bottom  profile  during  the  experiment  consisted  of  a  mild  offshore  slope  with  a  well 
developed  bar  approximately  100  m  from  shore  followed  by  a  trough  and  steeper 
foreshore.  (Figure  1).  Fluctuations  of  bottom  topography  in  the  surf  zone  were  minor 
throughout  the  data  collection  period  with  the  exception  of  the  passage  of  the  largest 
storm  that  moved  the  bar  80  meters  seaward. 


IH.  OBSERVED  SPECTRAL  ENERGY  BALANCE 

The  spectral  energy  balance  of  shoaling  and  breaking  waves  observed  for  a  wide 
range  of  conditions  is  illustrated  here  with  four  case  study  examples  (Figure  2),  each 
based  on  a  1-hour  long  data  record.  For  each  case,  the  evolution  of  the  spectrum  E(f) 
across  the  beach  is  shown  in  Figure  3,  and  a  representative  video  image  of  the  sea  surface 
in  the  instrumental  area  is  shown  in  Figure  4.  In  case  I,  small  amplitude  swell  propagated 
across  the  instrumented  transect  with  little  breaking  (Figure  4).  The  evolution  of  the 
narrow  spectrum  shows  the  expected  growth  of  harmonic  peaks.  Breaking  and  associated 
energy  losses  occurred  primarily  on  the  beach  face  shoreward  of  instrument  array  (Figure 
4). 

Case  studies  II  and  III  are  from  the  two  major  storm  that  passed  through  the 
region  during  the  experiment.  In  both  cases  the  irregular  surf  zone  with  spilling  breakers 
extended  across  the  entire  350m-long  transect  (Figure  4).  The  evolution  of  the  spectra  of 
these  broad  wind  seas  across  the  surf  zone  is  similar  with  a  strong,  nearly  uniform  decay 
of  spectral  levels  in  the  energetic  part  of  the  spectrum  (Figure  3).  In  case  IV,  during  the 
waning  of  the  second  Noreaster,  moderately  energetic  swell  was  observed.  This  case  is 
characterized  by  a  narrow  surf  zone  confined  to  the  inner  part  of  the  instrument  array, 
with  some  intermittent  breaking  occurring  farther  offshore  on  the  sand  bar  (Figure  4).  The 
spectra  observed  in  case  IV  indicate  a  pronounced  second  harmonic  peak  outside  the  surf 
zone.  Inside  the  surf  zone  both  the  primary  spectral  peak  and  the  second  harmonic  peak 
are  attenuated. 

To  examine  the  role  of  nonlinear  interactions  and  dissipation  processes  in  the 
observed  spectral  evolution  of  these  four  case  studies,  estimates  of  the  energy  flux 
gradient  Fx(f)  (Eq.  3)  and  the  nonlinear  transfer  SJf)  (Eq.  4)  are  compared  in  Figures  5-8. 
The  dissipation  term  S^fjin  the  energy  balance  can  be  inferred  as  the  difference  between 
the  estimated  Fx(f)  and  S^f)  (Eq.  1).  In  the  non-breaking  wave  Case  I  (Figure  5)  the 
observed  values  of  S^  and  Fx  are  approximately  equal  at  all  locations  except  for  the 
innermost  two,  indicating  that  the  observed  spectral  changes  (Figure  3)  are  the  result  of 
nonlinear  energy  exchanges,  consistent  with  Boussinesq  model  results  for  similar 


conditions  presented  in  Elgar  et  al.,  (1997)  and  Norheim  et  al,  (1998).  The  positive 
values  of  Sn,  and  Fx  in  the  range  0.1-0.4  Hz  indicate  a  transfer  of  energy  from  the 
dominant  0.07  Hz  swell  to  high-frequency  harmonic  components,  consistent  with  the 
Boussinesq  model  predictions  of  Elgar  et  al.  (1997)  and  Norheim  et  al,  (1998).  The 
largest  transfers  occur  slightly  inshore  of  the  bar  crest  (H,I)  and  on  the  beach  face  (L). 
Comparisons  of  Fx  and  S^  at  low  frequencies  (<0. 1  Hz)  are  not  shown  because  partial 
reflection  of  the  small  amplitude  swell  from  the  beach  face  (Elgar  et  al.,  1997),  caused 
large  errors  in  the  energy  flux  estimates  (i.e.,  standing  wave  nodes  and  antinodes). 
Reflection  effects  are  apparent  in  the  large  fluctuations  of  Fx  observed  at  the  innermost 
stations  K  and  L  (Figure  5). 

In  case  II  with  energetic  seas  breaking  across  the  entire  transect  (Figure  6),  the  S^ 
estimates  are  characterized  by  a  negative  lobe  centered  at  the  spectral  peak  frequency,  and 
a  broader  positive  lobe  at  higher  frequencies.  The  integral  of  S^f)  over  all  frequencies 
is  zero  because  the  triad  interactions  conserve  energy.  The  S^  estimates  show  that  energy 
is  primarily  transferred  from  the  energetic  part  of  the  spectrum  to  higher  frequencies.  The 
observed  energy  transfer  increases  from  the  outer  most  station  C  (A  and  B  were  not 
operational  during  this  event)  to  a  maximum  value  on  the  sand  bar  crest  (F)  decreasing 
to  small  values  near  the  shoreline  (i.e.,  where  most  of  the  energy  has  been  dissipated). 
The  estimated  energy  flux  gradients  Fx  show  a  negative  lobe  in  the  energetic  part  of  the 
spectrum  that  nearly  equals  the  negative  lobe  in  the  S^  estimates.  The  surprisingly  close 
balance  of  Fx  and  S^  demonstrates  that  the  large  energy  losses  observed  in  the  energetic 
part  of  the  spectrum  (Figure  3)  are  accounted  for  by  nonlinear  transfers  to  higher 
frequencies,  and  thus  no  significant  dissipation  takes  place  in  the  energetic  part  of  the 
spectrum. 

The  observed  energy  flux  gradients  and  nonlinear  transfers  do  not  balance  at 
higher  frequencies.  Whereas  SJf)  estimates  show  a  broad  positive  lobe  above  the  spectral 
peak  frequency,  the  associated  Fx(f)  estimates  are  close  to  zero.  These  results  suggest  that 
the  energy  transferred  through  nonlinear  interactions  to  higher  frequencies  cannot  be 
absorbed  in  the  high-frequency  tail  of  the  spectrum  and  is  presumably  dissipated, 


qualitatively  consistent  with  the  saturation  of  wave  spectra  in  the  surf  zone  observed  by 
Thornton  (1977,  1979).  In  the  most  energetic  Case  HI,  similar  fully  developed  seas  with 
a  lower  peak  frequency  (.09  Hz)  were  observed  (Figures  3,4).  In  this  severe  storm,  several 
instruments  failed  and  thus  the  energy  balance  could  not  be  evaluated  on  the  seaward  side 
of  the  bar  crest.  Nevertheless  the  estimates  of  Snl(f)and  Fx(f)  (Figure  7)  are  remarkably 
similar  to  the  case  II  estimates,  with  a  close  balance  between  energy  losses  in  the 
energetic  part  of  the  spectrum  (the  negative  lobe  of  Fx(f))and  energy  transfers  to  higher 
frequencies  (the  negative  lobe  of  S^f)),  and  an  apparent  saturation  of  the  spectrum  at 
high  frequencies  (S^  >  0  and  Fx  ~  0).  As  in  Case  II  the  largest  energy  transfers  are 
observed  at  station  F  on  the  bar  crest.  At  the  innermost  stations  K  and  L,  fluctuations  in 
Fx(f)  at  lower  frequencies  indicate  that  wave  reflection  from  the  steep  beach  face  (i.e., 
nodes  and  antinodes)  caused  significant  errors  in  the  energy  flux  estimates. 

A  few  days  later  (Case  IV)  the  storm  has  moved  offshore  and  swell  with  a  peak 
frequency  of  0.09  Hz  (Figure  3)  propagated  though  the  array.  Seaward  of  the  bar  crest 
at  stations  A  and  B  energy  is  transferred  from  the  spectral  peak  to  it's  secondary  harmonic 
and  higher  frequencies  (Figure  8).  The  negative  values  of  Fx  observed  at  station  B  do  not 
balance  the  positive  S^  estimates,  suggesting  some  energy  is  dissipated.  No  data  was 
collected  on  the  bar  crest  but  the  large  reduction  in  spectral  levels  between  stations  B  and 
F  (Figure  3)  indicates  strong  dissipation  takes  place  in  intermittent  wave  breaking  on  the 
bar  crest.  Between  the  bar  and  the  beach  face  at  stations  F-J,  nonlinear  energy  transfers 
are  approximately  balanced  by  the  estimated  energy  flux  gradients.  Whereas  at  station  F, 
energy  is  transferred  from  the  primary  peak  to  the  2nd  harmonic,  further  inshore  at  stations 
I  and  J  energy  is  transferred  from  the  second  harmonic  to  higher  frequencies.  Finally,  at 
the  shallowest  stations  K  and  L  located  well  inside  the  surf  zone,  the  S^  estimates  show 
large  transfers  of  energy  from  the  energetic  part  of  the  spectrum  to  higher  frequencies. 
As  in  cases  II  and  III,  the  Fx  estimates  closely  balance  S^  in  the  energetic  part  of  the 
spectrum,  and  Fx  ~  0  at  high  frequencies,  indicating  that  the  decay  of  the  energetic  part 
of  the  spectrum  results  from  nonlinear  transfers  to  the  saturated  high-frequency  tail  of  the 
spectrum  where  the  energy  is  dissipated. 
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IV.  DISCUSSION 

The  nonlinear  transfer  of  energy  to  or  from  wave  components  with  frequency  f  is 
controlled  by  the  phase  0(f)  of  the  bispectrum  integral  b(f)  (Eq.  2).  In  intermediate  water 
depths  (kh=0(l))  the  theoretical  value  of  this  so  called  biphase  is  0°  or  180°  (Hasselmann 
et  al.,  1963).  For  6  =  0  or  180°,  the  imaginary  part  of  b(f)  =  0  and  thus  no  transfers  of 
energy  take  place.  At  these  depths  the  bispectrum  describes  the  nonresonant  coupling 
between  the  primary  wind  waves  and  bound  secondary  waves.  In  shallow  water  (kh«l) 
biphases  evolve  allowing  for  energy  transfers  in  near-resonant  triad  interactions.  Negative 
biphase  values  correspond  to  a  loss  of  energy  for  components  with  frequency  f,  and 
positive  values  (O<0(f)<18O°)  indicate  energy  gains.  The  biphase  0  also  characterizes  the 
shape  of  waves  that  typically  evolve  from  symmetric  profiles  in  deep  water  (0=0  or  180°) 
to  pitched-forward  crests  prior  to  breaking  in  shallow  water  (0=  ±  90°)  (Elgar  and  Guza, 
1985). 

Biphase  estimates  are  shown  in  Figure  9  for  all  four  case  studies  at  locations: 
seaward  of  the  bar,  near  the  bar  crest,  and  on  the  beach  face.  Seaward  of  the  bar,  the 
biphases  are  close  to  the  intermediate  depth  values  0°  or  180°.  Case  I  shows  more  scatter 
probably  because  the  nonlinearity  is  weak  in  this  case  and  thus  the  bispectral  estimates 
have  greater  statistical  uncertainty.  Near  the  bar  crest,  the  bisphases  have  evolved  in 
cases  I-1TJ  to  about  135°  near  the  peak  frequency  and  to  30-90°  at  higher  frequencies.  The 
strongest  asymmetry  is  observed  in  the  non-breaking  Case  I  where  0  is  close  to  90°  at  the 
second  harmonic  peak  frequency.  Interestingly  in  cases  I-III  biphases  observed  on  the 
beach  face  show  a  relaxation  towards  the  values  0  and  180°  for  symmetric  wave  profiles. 
On  the  other  hand,  in  case  IV  biphases  observed  inshore  of  the  bar  crest  (where  waves 
reform  after  breaking  on  the  bar)  are  close  to  0  and  180°,  and  strong  asymmetry  does  not 
develop  until  waves  break  again  on  the  beach  face.  These  results  indicate  that 
asymmetries  in  wave  profiles  are  most  pronounced  during  or  before  the  onset  of  wave 
breaking  and  disappear  in  the  inner  region  of  wide  surf  zones. 

The  energy  balance  observed  in  case  studies  II-IV  with  energetic  breaking  waves 
(i.e.,  the  surf  zone  extends  seaward  of  the  sand  bar)  is  qualitatively  similar  with  nonlinear 
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interactions  transferring  energy  from  the  energetic  part  of  the  spectrum  to  higher 
frequencies  where  the  energy  is  dissipated.  This  surf  zone  balance  is  shown  in  Figures  10 
and  1 1  to  hold  for  all  observations  with  an  energetic  wide  surf  zone  during  the  two  month 
long  experiment.  Only  77  low-tide  observations  with  incident  wave  variances  measured 
in  8  m  depth  >  2000  cm2  (i.e.,  significant  wave  heights  >  1.8  m)  are  included  to  restrict 
the  comparisons  to  wide  surf  zones  that  extend  across  the  entire  instrumented  transect. 
Estimates  of  Fx(fp)  and  S^fp)  at  the  spectral  peak  frequency  (fp)  are  compared  at  locations 
A-L(midway  between  adjacent  instruments)  in  Figure  10. 

Between  the  seaward  end  of  the  transect  and  slightly  inshore  of  the  bar  crest 
(locations  A-I),  S^fp)  and  Fx(fp)  are  both  negative  and  approximately  equal.  Correlation 
coefficients  between  S^  and  FT  vary  between  .76-. 96  with  least-squares  fit  lines  that  are 
close  to  (within  5-30%)  a  one-to-one  correspondence.  These  estimates  demonstrate  that 
the  large  energy  losses  observed  in  the  surf  zone  at  the  spectral  peak  frequency  are 
accounted  for  by  nonlinear  transfers  to  higher  frequencies  rather  than  direct  dissipation 
of  wave  energy  in  the  energetic  part  of  the  spectrum. 

Gose  to  shore  (locations  J-L)  the  Snl(fp)  estimates  are  small  and  in  poor  agreement 
with  Fx(fp)  estimates.  At  these  inner  surf  zone  locations,  energy  levels  may  be  reduced 
enough  that  significant  reflected  components  from  the  beach  (e.g.,  an  antinode  near  the 
shoreline)  contribute  large  errors  to  the  energy  flux  gradient  Fx(fp)  estimates  (e.g.,  the 
large  positive  Fx(fp)  estimates  at  L). 

To  characterize  the  energy  balance  at  higher  frequencies,  the  S^f)  and  Fx(f) 
integrated  over  the  frequency  range  of  .3-. 5  Hz  are  compared  in  Figure  11.  Whereas  the 
S„i  estimates  show  consistently  positive  values  (i.e.,  energy  transfers  to  high  frequencies, 
Fx  estimates  are  comparatively  small.  These  results  confirm  that  the  high-frequency  tail 
of  the  spectrum  is  saturated  in  the  surf  zone  and  the  energy  transferred  to  high 
frequencies  in  nonlinear  interactions  is  immediately  dissipated.  The  largest  energy 
transfers  (and  hence  dissipation  rates)  are  observed  at  stations  F-L,  near  the  crest  of  the 
sand  bar,  and  L  on  the  beach  face  where  the  strongest  breaking  activity  is  expected. 


12 


V.  SUMMARY  AND  CONCLUSIONS 

The  spectral  energy  balance  of  shoaling  and  breaking  waves  was  examined  with 
two  months  of  field  observations  from  a  350m-long  cross-shore  transect  of  14  pressure 
sensors  deployed  on  a  barred  ocean  beach  near  Duck,  NC  as  part  of  the  Duck  94 
Experiment  (Elgar  et  ai,  1997).  The  energy  flux  gradient  term  Fx(f)  and  the  nonlinear 
source  term  S^f)  in  the  energy  balance  equation  were  estimated  at  13  cross  shore 
positions  from  approximately  1400  hour-long  pressure  records.  The  energy  flux  gradient 
Fx(f)  was  estimated  from  the  measured  changes  in  the  wave  spectrum  E(f)  between  two 
adjacent  instruments.  The  net  nonlinear  energy  transfer  SJf)  resulting  from  near-resonant 
triad  interactions  was  estimated  from  the  observed  wave  bispectrum  B(fl5f2)  based  on  a 
stochastic  Boussinesq  model  (Herbers  and  Burton,  1997). 

In  low-energy  conditions,  when  waves  propagate  across  the  instrumented  transect 
with  little  or  no  breaking,  the  observed  small  nonlinear  energy  transfers  to  higher 
frequencies  are  approximately  balanced  by  the  growth  of  the  spectrum  at  high  frequencies, 
consistent  with  earlier  studies  Elgar  et  <z/.(1997),  Norheim  et  a/.(1998).  In  high-energy 
conditions  with  waves  breaking  across  the  entire  transect,  the  S^f)  estimates  consistently 
show  a  large  negative  lobe  at  the  spectral  peak  frequency  and  a  broad  positive  lobe  at 
higher  frequencies,  consistent  with  a  strong  nonlinear  transfer  of  energy  from  the  spectral 
peak  to  a  broad  range  of  higher  frequencies.  The  largest  energy  exchanges  are  observed 
on  the  crest  of  the  sand  bar,  where  wave  breaking  is  most  intense.  The  observed  negative 
energy  flux  gradients  in  the  energetic  part  of  the  spectrum  closely  balance  the  nonlinear 
transfers.  Theses  observations  show  that  the  large  energy  losses  in  the  spectral  peak  of 
breaking  waves  are  accounted  for  by  nonlinear  transfers  to  higher  frequencies,  and  thus 
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no  significant  dissipation  occurs  in  the  energetic  part  of  the  spectrum.  At  higher 
frequencies  the  observed  Fx  are  small  and  do  not  balance  the  positive  S^  values.  These 
results  support  earlier  observations  (Thornton,  1977,  1979)  that  wave  spectra  in  the  surf 
zone  are  saturated  at  high  frequencies  (i.e.,  energy  transferred  through  nonlinear 
interactions  cannot  be  absorbed  in  the  high-frequency  tail  of  the  spectrum  and  is 
presumably  dissipated). 

In  conclusion,  the  present  observations  of  the  spectral  energy  balance  in  the  surf 
zone  show  that  near-resonant  triad  interactions  transfer  energy  from  the  energetic  part  of 
the  wave  spectrum  to  higher  frequencies  where  the  energy  is  dissipated. 
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APPENDIX 

Figure  1.  Instrument  locations  in  the  Duck94  experiment.  Squares  indicate  the 
14  cross-shore  positions  of  a  co-located  pressure  sensor  and  sonar  altimeter.  A 
(outermost)  through  M  (innermost)  indicate  the  midway  points  between  adjacent 
instruments. 

Figure  2.  Variability  of  incident  wave  variances  estimated  in  the  frequency  range 
0.05-0.25  Hz  from  bottom  pressure  measurements  in  8  meter  depth  (instruments  deployed 
and  maintained  by  the  U.S.  Army  Corps  of  Engineers)  using  a  linear  theory  depth 
correction.  Hourly  variance  estimates  are  shown  versus  time  for  the  two  months  of 
continuous  data  collected.  Dotted  vertical  lines  denote  the  four  case  studies  analyzed  in 
section  3. 

Figure  3.  Observed  shoaling  evolution  of  wave  spectra  in  the  four  case  studies. 
In  each  case  spectra  at  four  locations  spanning  the  transect  (Figure  1)  are  shown. 

Figure  4.  Representative  surf  zone  photographs  (courtesy  of  T.C.  Lippmann)  of 
four  case  studies  obtained  with  video  cameras  mounted  on  a  nearby  tower  about  50  m 
above  the  sea  surface.  The  top,  middle-left  and  lower-left  photographs  all  cover  the  same 
area  including  the  inner  half  of  the  instrumented  transect  (solid  line).  The  middle-right 
and  bottom-right  photos  provide  simultaneous,  slightly  overlapping  (~20m)  images  that 
cover  the  outer  half  of  the  instrumented  transect  (solid  line).  The  outward  looking  camera 
was  not  operational  during  cases  I  and  II. 

Figure  5.  Estimates  of  the  energy  flux  gradient  Fx(f)(dot-dash)  and  the  nonlinear 
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transferred  energy  S^f)  (solid)  versus  frequency  f  at  10  cross-shore  locations  in  case 
study  I.  The  beach  profile  and  the  cross-shore  locations  where  the  energy  balance  is 
evaluated  are  indicated  in  the  bottom  right  panel.  Arrows  indicate  the  spectral  peak 
frequency. 

Figure  6.  Estimates  of  Fx(f)  and  S^f)  in  Case  study  II  at  10  cross-shore  locations 
(same  format  as  figure  5). 

Figure  7.  Estimates  of  Fx(f)  and  Snl(f)  in  Case  study  III  at  8  cross-shore  locations 
(same  format  as  figure  5). 

Figure  8.  Estimates  of  Fx(f)  and  Snl(f)  in  Case  study  IV  at  9  cross-shore  locations 
(same  format  as  figure  5). 

Figure  9.  Biphase  0  versus  frequency  f  (Hz)  seaward  of  the  bar,  near  the  bar 
crest,  and  on  the  beach  face.  From  top  to  bottom:  case  I  at  locations  C,G,L,  case  II  at 
C,G,L,  case  III  at  A,F,L,  and  case  IV  at  A,F,L,. 

Figure  10.  Nonlinear  energy  transfer  S^  versus  the  energy  flux  gradient  Fx  at  the 
spectral  peak  frequency  in  energetic  surf  zone  conditions.  Each  scatter  diagram  shows  the 
comparison  at  one  cross-shore  location  (comparisons  at  locations  D  and  E  where  only  few 
observations  were  collected  are  not  shown).  Each  dot  represents  an  estimate  based  on  a 
1 -hour-long  data  record.  A  linear  regression  line  (dash-dot)  and  correlation  coefficient 
between  Fx  and  S^  are  indicated  in  each  scatter  diagram.  For  reference,  a  representative 
beach  profile  with  the  instrument  locations  is  shown  in  the  lower  right  panel. 

Figure  11.  Nonlinear  energy  transfer  S^  versus  the  energy  flux  gradient  Fx,  both 
integrated  over  the  high-frequency  range  0.3-0.5  Hz,  in  energetic  surf  zone  conditions 
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(same  format  as  figure  10). 
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